Analysis of path integrals at low temperature: 
Box formula, occupation time and ergodic approximation 
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We study the low temperature behaviour of path integrals for a simple one-dimensional model. 
Starting from the Feynman-Kac formula, we derive a new functional representation of the density 
matrix at finite temperature, in terms of the occupation times of Brownian motions constrained to 
stay within boxes with finite sizes. From that representation, we infer a kind of ergodic approx- 
imation, which only involves double ordinary integrals. As shown by its applications to different 
confining potentials, the ergodic approximation turns out to be quite efficient, especially in the 
low-temperature regime where other usual approximations fail. 

■ PACS numbers: 05.30.-d Quantum statistical mechanics; 05.40.Jc Brownian motion 
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The knowledge of the density matrix at finite temperature T for few interacting particles, is of paramount impor- 
tance for studying equilibrium properties of quantum many-body systems. Indeed, such density matrices occur in 
low-fugacity expansions, while they are also key ingredients in quantum Monte Carlo simulations. In this context, 
the Feynman-Kac (FK) representation of the density matrix in terms of path integrals P, Q is particularly useful. 
On the one hand, it has been used for deriving exact (analytical) expressions for simple models @- On the other 
hand, beyond the well-known Wigner-Kirkwood expansion El around the classical limit, various approximations, 
non-perturbative in fi, have been introduced within that framework: for instance, the celebrated semi-classical ap- 
proximation (widely used in numerous situations), the variational approach of Feynman and Kleinert 0, or the 
^-H ' renormalized Wigner-Kirkwood scheme. 

^ ' TT, , , , 

The FK representation is described in Sec. [0] For the sake of simplicity, we consider the simple case of a single 
particle in one dimension submitted to a stationary confining potential. The mean spatial extension of the brownian 
paths, which intervene in the FK representation, is controlled by the thermal de Broglie wavelength of the quantum 
particle. At high temperatures, since that extension vanishes, exact asymptotic expansions, as well as direct numerical 
estimations, become rather straightforward. At the analytical level, Wigner-Kirkwood fi 2 -expansions around the 
classical limit, which have been derived long ago in other frameworks, are easily recovered. On the computational 
side, functional integration over paths can be replaced, within a good accuracy, by ordinary integrals over a discrete 
set of parameters, via the use of either free propagators in Trotter formula [g or simple modelizations of brownian 
paths 9]. However, at low temperatures, similar calculations loose accuracy, while they also become rather time- 
consuming. Moreover, the corresponding asymptotic structure of the density matrix, provided by the groundstatc 
contribution in its spectral expression, does not clearly emerge from its FK representation. In fact, in that temperature 
regime, the average extension of paths diverge, whereas the main contribution to the functional integral arises from 
paths with a finite extension, of the order of the localization length of the ground state: contributions of such minority 
paths are significantly different from the average contribution. Most of our knowledge is negative, i.e. tells us which 
trajectories are not important [Tol ]. so, in some sense, we will try to have a positive attitude 

In this paper, we derive a new functional representation of the density matrix, which is more suitable than the genuine 
FK formula for tackling the low-temperature regime (see Sec. Illlf) . The starting central observations are described 
in Sec. III. A. First, in the FK functional integral, only marginal paths with finite extensions (i.e. large deviations 
with respect to the average) contribute when T goes to zero. Second, many of such paths with quite different jagged 
shapes, provide similar contributions, mainly determined by the corresponding local occupation times. Thus, it is 
quite natural to collect paths into sets defined by their spatial extension and their local occupation times. That 
procedure allows us to transform exactly FK representation <|10|) into the so-called box formula l|35|l . That formula 
is defined via the introduction of paths constrained to stay in a box with size i and characterized by an intermediate 
flight time s (in (3h units with (5 = l/fc^T). It involves a double ordinary integral over I and s, combined to functional 
integrals over local occupation times associated with the constrained paths. As required, box formula provides 
a better understanding of the low-temperature behaviour of the density matrix than FK representation i|10|) . When 
T vanishes, leading contributions obviously arise from typical sizes I much smaller than de Broglie wavelength. In a 
forthcoming paper, we will argue how groundstate quantities emerge from box formula l|35|) . by using scaling properties 
of the probability distribution function (PDF) of occupation times at low temperatures. 
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Beyond its conceptual interest for understanding low-temperature behaviours of path integrals, box formula 1351) 
also allows us to derive new approximations. This is illustrated in Section llVI where we present the so-called ergodic 
approximation. That approximation results from the truncation to first order of cumulant expansions of the functional 
averages. It amounts to replace each local occupation time associated with a given constrained path, by its average over 
the PDF of such occupation times. This can be viewed as some kind of ergodic hypothesis, because the (imaginary) 
time average of the potential experienced by the particle along that path, is then replaced by a spatial average with 
a measure defined by the previous mean occupation-time. Ergodic expression l|39|) involves only an ordinary double 
integral over £ and s. The key ingredient, namely the mean occupation-time (in units of l/l), depends on three 
dimensionless parameters. Using its low and high temperature behaviours (derived analytically), we propose simple 
tractable expressions for that quantity, which turn out to be quite accurate at any temperature. 

Section[V]is devoted to the applications of the ergodic approximation to various simple forms of the confining potential. 
First, we determine the asymptotic analytical forms of the approximate density matrices at both low and high 
temperatures. When T diverges, Wigner-Kirkwood expansion is partially recovered. When T vanishes, the main 
features of the exact behaviours are well reproduced. The approximate density matrices then do factorize as a product 
of a Boltzmann factor associated with a given energy, times a function of position only: this provides satisfatory 
approximate expressions for the groundstate energy and wavefunction. Second, numerical calculations are performed 
at finite intermediate temperatures. As expected from the previous analytical results, the ergodic approximation 
turns out to be quite reasonably accurate over the whole range of considered temperatures, and discrepancies with 
(numerically) exact results (in part inferred from the spectral representation) do not exceed a few percent. Moreover, 
it significantly improves over the well-known semi-classical approximation, which completely fails at low temperatures 
(except of course for the harmonic potential where semi-classical calculations become exact). 

Further applications and extensions of our approach will be presented elsewhere. In particular, the ergodic approxi- 
mation can be applied to unconfining potentials, and it can be straightforwardly extended to higher dimensions and 
to systems with several particles. Furthermore, other approximations might be easily introduced by starting from box 
formula as argued at the end of Section III. 



II. PATH INTEGRAL FRAMEWORK 



In this Section, we first define the model and then introduce the FK representation of the corresponding density 
matrix. Next, we briefly recall the efficiency of FK formula for describing the high-temperature regime, and we argue 
about its drawbacks for analyzing low-temperature behaviours. 



A. The model 



We consider a quantum non-relativistic particle of mass m in one dimension z, submitted to a confining potential V, 
i.e V{z) — > oo when \z\ — > oo. Its Hamiltonian reads 

h 2 d 2 
2m az z 

Introducing inverse temperature /3 — l/(ksT) with Boltzmann constant hs, we define the Gibbs operator 

p = e X p[-/3H], (2) 
the matrix elements of which define the so-called density matrix p(x, y, (3) = (x \p\ y). The partition function 

/+oo 
dzp(z,z,/3) (3) 

is well behaved thanks to the confining nature of potential V (in other words, integral over z does converge because 
p(z,z,{3) decays sufficiently fast at large distances). The (normalized) probability density to find the particle at 
position x then reduces to 

if, l R\ P( X > X >P) (A\ 
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For further purposes, it is also convenient to introduce the spectral representations of density matrix and partition 
function, which read 

+00 

p(x,y,(3) =J2Mx)<t>t(y)exp(-pE k ) (5) 

and 

+00 

Z(0)=Y,™p(-(3Ek). (6) 

k=0 

In such representations, Ek is the fc-th eigenvalue of H, while <f>k is its associated eigenfunction which satisfies 
Schrodinger equation 

i~^ + V{z)]Mz)=E k Mz). (7) 



B. Feynman-Kac formula 

Path integrals were first introduced for representing the matrix elements of the evolution operator exp(—iHt/h) 
associated with Schrodinger equation It was soon realized that a similar path integral representation for matrix 
elements of Gibbs operator at inverse temperature (3 can be inferred via the formal substitution t — > —i[3h. The 
latter representation for p{x,y,(3) involves an integral over all paths uj(-) going from lo(0) = x to uj{[3K) = y in 
a time (3h. Moreover, the corresponding integrand is the exponential factor exp[— S(lv(-))/H], where S(lu(-)) is the 
action associated with path ui(-). That path integral is proved to be mathematically well defined for a large class 
of potentials. Its convergence is indeed ensured, roughly speaking, by the Gaussian-like decay of the kinetic part of 
exp[— S(cu(-))/h] for large paths (on the contrary, path integrals for exp(~iHt/h) arc, in general, ill-defined because 
they involve sums of oscillating phase factors). 

The genuine path integral representation of p(x, y, (3) can be rewritten within the parametrization lu(u) = (1— s)x+sy+ 
Ad£(s), where s — u/((3K) is the dimcnsionlcss time in (3h units, while £(s) is a Brownian bridge satisfying boundary 
conditions £(0) = £(1) = and Xd = (f3h 2 /to) 1 / 2 is the de Broglie wavelength. On the one hand, the kinetic part 
of exp[— S(uj(-))/h] then provides normalized Wiener measure 2?w(£) associated with the Brownian bridge process. 
That Gaussian measure is entirely defined by its first two moments 

! m = 0, (8) 

Jn 

/ 2V(0 =min(si,s 2 )(l-max(si,S2)). (9) 

J a 

On the other hand, the potential part of exp[— S(iv(-))/h] reduces to the Boltzmann factor associated with time average 
of V along path (1 — s)x + sy + \d£(s). The resulting so-called Feynman-Kac formula reads |2| 

p(x, y, (3) = eXph(x Z- y)2/(2A ^ )] / 2V(0 exp (-(3 [ ds V ((1 - s)x + sy + Xd^s))) (10) 
V Z7r \ D Jn \ Jo J 

where f2 = {£(•)} is the infinite set of realizations of the Brownian bridge process. 

In the Feynman-Kac formula (| 1 - the Wiener measure is intrinsic to Brownian motion and does not depend on any 
physical parameter. Expression Q of the covariance implies that the typical extension of a Brownian bridge is always 
of order 1. Mass to of the particle, as well as Planck's constant H only intervene in the de Broglie wavelength Ap. 
That length controls the size of quantum fluctuations, as shown by specifying l|10|) to diagonal elements, i.e. 

p(x, x, !3) - f V W (0 exp (-f3 [ dsV(x + X D £( S ))) . (11) 

V2ttXd Jn \ Jo J 

In time-average J* ds V (x + Ad^(s)), particle experiences the potential around position a; on a length scale obviously 
determined by Ad, as illustrated in Figure ^ 
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FIG. 1: Typical path which starts and ends at position x, along which particle experiences potential V. Dashed lines indicate 
the two edges of a fictitious box determined by the extremal deviations £1 and £2 of the path. Typical sizes of £1 and £2 are 
controlled by Ad 



Besides its interest for analytical or numerical calculations (see below), Feynman-Kac representation is also quite 
useful for mapping a quantum system into a classical one [Sill El- Indeed, formula (JTTJ can be reinterpreted in 
terms of a classical extended object, often called loop or polymer, denned by its position x and its internal shape 
parametrized by Xd£,(-)- Then, the r.h.s. of is nothing but the statistical average of a Boltzmann factor for the 
loop over its all possible shapes distributed according to Wiener measure. That interpretation, which holds for any 
quantum system with an arbitrary number of particles, is the key starting point of above mapping. 



C. High and finite temperatures 



At high temperatures, the de Broglie wavelength is small, so only paths which remain close to reference position x 
contribute significantly to p(x, x, (3) in (|TT1) . Therefore, the time-average of V along path x + XdS, (s) can be performed 
by replacing V {x + A_d£(s)) by its Taylor series around V (x) 



V(x + \ D Z(s))=V(x) + J2 



(x D g s )rd n v 

n\ dx n 



(12) 



(we assume that V is infinitely diffcrentiable). If we keep only the first term in (|12|l . which does not depend on £(s), 
we immediately recover the familiar classical form of p(x, x, [3) 



■exp[-/3V(i)], 



(13) 



where we have used normalization condition J^T^wiS.) = !■ Since Xd is proportional to h, quantum corrections to 
(f 1 31) are easily generated by a (formal) cumulant expansion of contributions associated with terms n > 1 in (|12ll . The 
resulting moments of Wiener measure are then readily computed by Wick's theorem which applies thanks to Gaussian 
character of T>w(Q- This allows us to retrieve the well-known Wigner-Kirkwood expansion in powers of h 2 



p{x,x,p) 



-0V{x) 

2vrAn 



12to da; 2 24m \dx J 



(14) 



Expansion (|14fl involves only even powers of H, because all odd moments of gaussian Wiener measure vanish by virtue 
of JHJ). That expansion is expected to be only asymptotic. Present derivation within FK representation turns out 
to be much easier than usual approaches within Wigner-distribution formalism (see e.g ref. |fj). It might be used 
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for calculating higher-order terms beyond ft 2 -correction (terms up to order h 6 have been already determined in the 
literature [lffi. Moreover, it clearly emphasizes that Q14[l is appropriate when the de Broglie wavelength is much 
smaller than the typical length of variation of the potential. Then, truncation of Ijl4|l up to a few terms should 
provide accurate values of the density matrix, as checked below. 

When the temperature decreases, the de Broglie wavelength increases, so Wigner-Kirkwood expansion is no longer 
applicable. At the analytical level, exact evaluations of the corresponding path integrals are not accessible in general, 
except for some simple models ,3]. Non-perturbative effects in h can be partially accounted for through various 
methods, among which the celebrated semi-classical approximation described further. The reliability of such methods 
remains questionable since they do not involve any small control-parameter like 114(1 . Nonetheless, for numerical 
purposes, their accuracy may be satisfactory, at least at not too low temperatures. 

If analytic informations on the density matrix become rather difficult to extract from FK representation, accurate 
values can be still readily obtained at finite temperatures via direct numerical estimations of FK path integrals. A 
first possible route starts with the (exact) convolution relation 

r N N 

p(x,x,0) = / Y[dz z Y[p(z t ,z l+ll l3/(N+l)) (15) 

•> i—l i—0 

with zq = z/v+i = x. For a sufficiently large value of N, each matrix element p(zi, Zj+i, /?/ (N + 1)) can be evaluated 
within FK formula (f 1 C)|> by neglecting the corresponding Brownian deviation A_d£(s) / (N + l) 1 / 2 in the argument of V. 
The resulting time average L dsV ((1 — s)x + sy) can be itself replaced, with a similar accuracy, by V((x + y)/2) 
(middle-point convention) for instance. The resulting approximate expression 

Jfl exp[ _ mN + 1} ± ({Zi + 2m)/2)] 

turns out to be quite accurate for rather small values of N, provided that Ad/(N + l) 1 / 2 is smaller than the typical 
length of variation of the potential. We stress that l|16f) becomes exact when N — > oo. In fact, continuous path 
integration with Wiener measure may be precisely defined in the framework of that limit process. 

A second route for direct numerical computations relies on the introduction of a set of N simple functions e n (s) 
(n = 1...N) with e n (0) = e„(l) = 0, such that they define a basis for Brownian bridges £(s) in the limit N — > oo 
|9| . Numerical evaluations of (| 1 II) are then carried out by replacing the functional integration over paths by ordinary 
integrals over coefficients o„ arising in the decomposition £(s) = a n e n (s). That method provides satisfactory results 
by using only a few number N of basis functions Similarly to the previous convolution method, it becomes exact 
when N — > oo. 



D. Low temperatures 

When temperature goes to zero, Ad diverges, so most paths explore a rather large region not restricted to the 
neighbourhood of reference position x. A direct analysis of FK expression (|11|) then becomes quite cumbersome. In 
particular, p(x, x, (3) must behave as 

p(x,x,(3) ^°° |0 o (x)| 2 exp[-/3£;o], (17) 

which immediately follows from spectral representation of density matrix ©. Indeed, the groundstate obviously 
provides the leading contribution in the r.h.s. of |JSj). The factorization of position and temperature dependences 
in (|17fl does not come out easily from FK formula l|llf> . where both position and temperature are coupled in an 
absolutely non-trivial way. By the way, notice that spectral representation JSJ is not well suited at intermediate or 
high temperatures, because one has then to sum a large number of contributions associated with various eigenstates 
(in particular, classical behaviour (|13f) cannot be straightforwardly recovered by using (JSJ). 

When temperature becomes too small, previous computational methods become more difficult to handle accurately. 
Even if some algorithmic tricks 16] may significantly reduce calculation-times within a satisfactory accuracy, they do 
not provide a better understanding of low-temperature behaviours. 
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III. BOX FORMULA 
A. The central observations 

Let us consider a symmetric confining potential V(z) with a minimum located at z = 0. On the one hand, low- 
temperature behaviour (|17H of the corresponding density matrix is mainly determined by the local shape of V{z) over 
finite length scale ao, which characterizes the spatial extension of the groundstate wave function <fio(x). On the other 
hand, for typical paths with size of order Ad, time-average potential f Q dsV (x + Ad£(s)) becomes, roughly speaking, 

of order J^ D dzV (z) /Ad when Ad is sufficiently large. At low temperatures, their contribution to the r.h.s. of l(TT|) 
is then controlled by the large- distance behaviour of V(z). Thus, low-temperature behaviour (fTTjl is not provided 
by typical paths. That argument can be implemented through semi-quantitative estimations for specific cases. If 



V{z) diverges as \z 
as exp [-c(/?) 1+n / 2 " 
exp[-j3J5 ]. 



(n > 0), contributions of typical paths to behave (discarding multiplicative powers of 0) 
with some positive constant c: they are exponentially smaller than leading Boltzmann factor 



The previous analysis suggests that, at low temperatures, leading contributions to the r.h.s. of Ijllfl arise from paths 
with a spatial extension of order ao, i.e. from quite small Brownian bridges with size |£(s)| of order oo/Ad- For 
such paths, time-average potential J Q ds V (x + Ad£(s)) is of order V{ciq), so the corresponding Boltzmann factor 
indeed is of order exp [—/3Eo]. Notice that, paths with very different shapes give raise to similar contributions, since 
the Wiener weights of the associated Brownian bridges remain of order (roughly speaking) exp(— J ds(£(s)) 2 /2) ~ 
exp(— ag/A|j) ~ 1. Consequently, when the temperature decreases, only a very tiny subset of paths gives a relevant 
contribution to the r.h.s. of (|1 1|) . In other words, important paths are not any more typical but, on the contrary, they 
can be viewed as large deviations. This explains why direct numerical evaluations of (|ll(l become rather difficult: the 
subset of important paths remains, in some sense, hidden in the entire configurational space. 

In order to extract from (|ll|l the relevant contributions at low temperatures, it is tempting to collect all paths with 
the same finite spatial extension, and then to sum over all possible extensions. That procedure is not easy to carry 
out directly in the r.h.s. of (|11|) . within a suitable partition of functional integration space over Brownian bridges. 
As described below, it is more convenient to transform first the density matrix within the operator representation, by 
introducing an auxiliary Hamiltonian which confines the particle inside a box with size £. 

B. The auxiliary Hamiltonian approach 

Let us introduce the auxiliary Hamiltonian 

H e = H° + V + V t , (18) 

where H° denotes the kinetic Hamiltonian. The additional potential Vg is defined by 

V e (z) = V [l-@(z + £) + e(z-l)}, (19) 

where O is the Heavyside function, while Vq denotes barrier height (Vb > 0). That potential tends to confine the 
particle inside a box with extension I. Expression (|18fl reduces to Hamiltonian when I goes to infinity. In that 
limit, for x kept fixed, diagonal part (x\pe \x) of density matrix pe = exp [—j3H(\ goes to (x|/3|x). 

The identity, 

/>OC 

{x\p\x) = (x\ PL \x) + dl {x\d e [p e ] \x) (20) 

which is valid for any reference extension L and for any height Vo of the confining potential is quite useful for our 
purpose. We set L = \x\, and we consider an infinitely high potential barrier Vb- Under that limit, the contact density 
goes to zero, 

lim (x\p {xl \x)=0, (21) 
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because the probability for the particle to stay on the boundary vanishes for infinitely high walls. Thus, identity 1201) 
is rewritten as 



(x\p\x) 



lim 

Vo — >oo 



d£ (x\dt \p t ] \x) 



|,-| 



The right hand side can be transformed by applying Dyson formula 

d £ \e- A W] = - f 1 d S e-^ A ^d e [A(£)} e~ sA ^\ 
L - 1 Jo 

valid for any operatorA(i?) which depends on I. This provides 



di [pi 



-(3 



ds 



— oo 



dz e-"* 1 — \z)(z\ d t [H t ] e-P* H ' 
dz e -W-')Bi \z) [S(z -£) + S(z + £)} (z\e-P' H * 



(3V f ds(e-^- s ^\l)(i\e-^ + e -^- s ^\ -£)(-£\e~^ Hi ) 



(22) 

(23) 

(24) 
(25) 
(26) 



from which we infer 

p(x,x,(3) = 



lim (3V 

Vq^OC 



At 



ds [ Pt (x, t, (3[l -s])pe (£, x, (3s) + p t (x, -£, (3[1 - s}) p t (-£, x, /3s)} . (27) 



That expression must be considered with some caution because when Vq diverges the different integrands vanishes. 
Therefore, additional work is necessary to control limit Vo — > oo. 

For that purpose, we rewrite the integrands in terms of constrained density matrix p® of the free particle submitted 
to confining potential V(_. According to the obvious equality 



/ m o , o\ Pt (xi,x f ,/3) 
Pi(x l ,xf,(3) = p i (x l ,x f ,(3)— — 

Pi \Xi, Xf, p) 



and since pe/p® remains finite when Vq — > oo, we define 



7 ± (x,£,s,/?) = lim (27T) 1 / 2 \ D ((3Vo)p° i (x,±£,f3s)p e (±£,x,l3(l-s)). 

Vo— >oo 



(28) 



(29) 



After using Feynman-Kac formula for both pi and p\ in ratio pe/pg, we transform l|27|l into 



p(x,x,(3) 



k! 

+ 00 



d£ / ds 



o 



2tt\ d 
9 + (x,£,s,P) 



f3s I du V{z{u)) 



exp 



(3s / duV (z{u)) 



n: 



exp 



exp 



/3(1 -s) du V{z{u)) 



(3(1 - s) / du V(z(u)) 



(30) 



Notation ( . ) u denotes an average over a constrained Brownian process which belongs to a set u. Paths z[u) are 
expressed in terms of Brownian bridges according to z(u) — x{l — u) ±u£ + ^/s\d£,(u) for . Moreover those paths 
must stay inside the box [—£,+£}, so that constraint defines the corresponding sets Clf. The statistical weight of a 
path in such a constrained average is its associated Wiener measure T>w{C)- For the sake of notational convenience, 
we do not explicitly write the dependences on both x and £ of Slf . 

The physical interpretation of functions g± clearly emerges from (I30f) . if we specify that general formula to the 
particular case V(z) = 0. Such functions are the (normalized) statistical weights of the constrained sets i7 ± — 
Q±\J£lf_ s , i.e. the' sum of statistical weights of all paths touching the boundaries of the box at time s (see Figure 
|2J). Therefore, in the following, we set g ± (x,£, s, (3) = g(J7 ± ). Analytical calculations of those weights are performed 
in Appendix A, by using the spectral representation of p®. 

The first step of our rewriting of Feynman-Kac representation, is achieved through formula l|30|l . Paths are indeed 
collected together according to their extension ±£. Notice that the touching time, s, is also crucial for defining the 
corresponding proper partition of the genuine integration space over all unconstrained paths. 




u=0 



u=l 



FIG. 2: Typical paths which start and end at position x. Solid and dashed lines represent paths which belong to fi + and to 
Q~ respectively. 



C. Introduction of averages over occupation times 



In a second step, we introduce the so-called occupation time, defined for each given path z(u) in £lf by 

9 z {x') = [ duS(x' -z(u)). 
Jo 



(31) 



The quantity 9 z (x')dx' is the total time passed in a neighborhood dx' of position x' by the particle when it follows 
Brownian path z(u). Of course, the total time passed in the whole box is always equal to 1, i.e. 



dx'9 z (x') = 1. 



(32) 



Time- averaged potential along path z(u) is then expressed in terms of occupation time via the obvious identity 

r+t 

duV(z{u))= dx'6 z (x')V(x'), (33) 

JO J-i 

valid for any path z{u). 

According to identity l|33|l . Boltzmann factors involved in averages 



exp 



-(3s I du V (z(u)) 



(34) 



only depend on the occupation time 8 z (x'). As illustrated in Fig. [31 various different paths may provide the same 
occupation time. As quoted in Section IIII Al their statistical weights are close together, so their contributions to 
average (13411 are almost identical. Thus, it is now tempting to collect all paths which provide the same occupation 
time 8(x'), via the introduction of the corresponding density measure T> n ±[8]. After expressing averages l|34|l over 



9 



Az(s) 




u=0 



u=l 



FIG. 3: Three different paths which start at x and end at +£, while the corresponding occupation times are identical. 



constrained paths, as averages over occupation times with (normalized) measure Z> n ±[0], we eventually obtain the 
box formula 



p(x,x,j3) = / dl/ds ^ IV n - [9] exp 

J\x\ JO 



'2-kXd 
>+ 



di [ ds ^ [ V n + [9] exp -/3s f dx'9(x')V (x') 
x \ Jo V27tAdJ s J-e 



/3s / dx'9(x')V(x') 
-i 

+i 



T> n - [9] exp 



(3(1 -s) dx'6{x')V{x') 



T> n + [9] exp 



+f 



0(1 -s) dx'9(x')V(x') 



.(35) 



If statistical weights g(Q ) are analytically known (see Appendix A), explicit expressions for density measures T> n ± [9] 
are not available. Such probability densities result from the summation of Wiener measures over all Brownian paths 
inside £lf which provide the same occupation time. That procedure is quite difficult to handle in closed analytical 
forms, and only the moments of T> n ± [9] can be computed explicitly. Nevertheless, box formula i|35|) turns out to be 
quite useful for various purposes, as suggested by the following simple comments and arguments. 

Contrary to the case of the genuine Feynman-Kac representation, leading contributions at low temperature merely 
emerge from box formula l|35l) . Indeed, for x of order ao, boxes with size £ of order a few ao do provide contributions of 
order exp(— (3E ), because weight factors g(fl^) are of order exp(—f3h 2 /(2m£ 2 )) for f<An (see Appendix A), while 

products of averages of Boltzmann factors over occupation times are of order exp(— (3 J Q dx'V (x 1 ) /£). That rough 
analysis will be implemented in a forthcoming paper, where we show more precisely how low temperature behaviour 
(|17fl arises from scaling properties of distributions T> n ± [9] in the regime £ <§; . 

If Brownian paths are quite noisy, the corresponding occupation times may have rather regular shapes. This is 
illustrated in Fig. which shows a regular path on the one hand, an a very jagged one on the other hand: both 
paths provide occupation times with regular shapes displayed in Fig. 0J> (notice that the corresponding contributions 
to Boltzmann factors are of course different). That observation leads to a first type of approximations based on 
simple modelizations of distributions T>^± [9] within restricted sets of elementary functions which represent the various 
occupation times (they will be described elsewhere). A second type of approximations is based on the truncation of 
cumulant expansions, where key ingredients do exhibit regular behaviours with respect to spatial positions. One of 
them is presented further in Section llVI 



Eventually, notice that box formula l|35|) can be immediately extended to any potential V(x), not necessarily confining. 
Similar formula can be derived in higher dimensions or for systems with more than one particle. Those various 
extensions may include different box shapes or different confining potentials Vq(x). 
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b) 



u=0 



X'l 






e- z (x'> \ 







FIG. 4: Panel (a) presents two different paths which start at position x at time u = and end at £ for time u = 1. Panel (b) 
shows their corresponding occupation times 6 z (x'). 



IV. ERGODIC APPROXIMATION 

A. Truncation of the cumulant expansion 

Formally, averages over distributions T> n ± [9] involved in the r.h.s. of box formula (|35fl can be represented by their 
infinite cumulant expansions. A natural approximation consists in truncating that expansion up to its first term, i.e. 
we replace the averages 



J V n ±[9]ex P 



j3s dx'9{x')V{x') 



by 



exp ^-/3s J dx' < 6{x') > n ± V (V)J 



(36) 



(37) 



We call that lowest order approximation ergodic, since it would be exact if paths were experiencing all parts of the 
potential with a probability independent of time. It amounts to replace (imaginary) time averages of the potential by 
spatial averages with a measure defined by mean occupation-time 



< 9(x') >, 



(38) 



Inserting Ij37|) into Ij35|) . we obtain the subsequent ergodic approximation for density matrix, 



d7/d S ^ex P 
x \ Jo \2-k\d 

+ d7/d S 4eiex P 
\ x \ Jo V27tA_d 



[3j^Ax'(s< 6{x') > n - +[1 - a] < 9(x') > n -_ J V (x') 
(3 p dx' (s < e(x') > n + +[1 - s] < 9{x') > Q +_)V{J) 



(39) 



B. Mean occupation time in terms of Brownian Green functions 



The average occupation time can be obtained within two different methods: either by using operator algebra (see 
Appendix B) or, in a more physical way, by using Brownian motion properties, as shown below. First we compute 
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the mean occupation time around position x', over Brownian paths at inverse temperature f3' which start from xi at 
time u — and end &t Xf for time u = 1. Those Brownian paths are constrained within the interval [—£, +£], and 
they define a set oj = {z(u) = Xi(l — u) + x fU + A£(it)} with A 2 = (3'ti 2 jm. 

Introducing the constrained probability density f^z^u) to find the Brownian particle at the position z for time u, 
the mean occupation time reads 



dufu,(x',u). 



(40) 



The probability density fu{x', u) can be generated by using the Green function G u (z, t\zo, to) of the diffusion equation, 
solved with Dirichlct boundary conditions G UJ (±£, t\zo, to) — and initial condition G ul (z, to\zo, to) = S(z — zq). That 
Green function reduces to 



G u (z,t\z ,t ) = J^2^n(z /^n(z/e) 



exp 



n=l 



(41) 



with ipn(z/ '£) = sin [wk(z + £)/(2£)]. Since Brownian motion is a Markov process, the constrained probability density 
fuj{x', u) is then 



, Guj{x', u\xj,0) Gu{xf, 1 - u\x', 0) 

U[X > U) - G u (x f ,l\xi,0) 



(42) 



Now, we focus on the specific case u> = Q s (i.e. Xf 
Gq- (x 1 , u\x, 0) can be rewritten as 



G Q -(x',u\x,o) = \m (7'7'^y^' m 1 ■ 



Xi = x and (3' = (3s). According to l|41|) . the Green function 

(43) 



where 771 is a dimensionless function depending on four dimcnsionless variables. Similarly, we find 



G Q -(-£,l-u\x',0) 
G n -(-£,l\x,0) 

If we introduce the dimensionless function $ defined as 



x x 1 Au-y/i 



(44) 



we obtain 



$(a',a,y) = / du 771(0:, a', y, u)rj 2 (a, a', y, u), 
Jo 



c' Xd\/s 



(«>nr - 7*17.7. 



(45) 



(46) 



Furthermore, thanks to the symmetry properties of Green functions, G llJ (z,t\zo,to) — G UJ (zo, t\z, to) = G u [— z,t\ — 
z ,t ), (6(x')) n + is also given by 



(d(x>)) 



—x —x' \d\/~S 
T' ~T' £ 



(47) 



Inserting those expressions of mean occupation times for various sets flf into l|39(l , the ergodic form of the density 
matrix is rewritten as 

p erg (x, x, p) — / d£ / ds — -= exp I — (3 I aa 

J\x\ Jo V27tA_d V J-i 

«+oo pi _f „ i) _ a\ / /•! 



2/' 



2f 



( , ; dg fl( ^5 s '^ exp(-/3 / da' 



2ttA 



D 



S <I> I — a' nXD ^~ S 



V (a'*) 



with (7(2;, s, /3) = g(flj). In l|48l) . all quantities are explicitly known in terms of simple scries involving Gaussian and 
trigonometric functions. Box weight g(x, £, s, (3) is computed in Appendix A, while a similar expression is derived for 
$ in the next subsection (see formula (|52|l 'l. Within the ergodic approximation, we are left with the evaluation of two 
ordinary integrals, which is of course much easier than a direct evaluation of the genuine functional integrals. 
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C. Analytical estimations for $ at low and high temperatures 

1. General expression for $ 

First, we derive a general expression for <E>, valid at any temperature. The spectral expressions of r\\ and r\i are 



rji(a,a',y,u) = ^ ip n (a)ijj n (a') exp 

n=l 
oo 

V^n (a')n exp 



2/ 2 2 



m (a,a',y,u) 



v 2 

— n 1 — it 

2 1 J 



n exp 



V 2 2 

n 

2 



(49) 



(50) 



Inserting those expressions into formula (|45fl . we find the general expression for 

oo 



$(a,a',y) 



y 2 2 

n 

2 



2 00 exp 
"a V'n(a)V'r 1 (a')' ( /'fe(a') fc 



exp 



» 1,2 



fc 2 — n 2 



X>*( 



ajnexp 



n=l 



?/ 2 2 

— n 
2 



(51) 



Normalization condition Jjt $(q!, a', y) = 1, which follows from l|32(l . is indeed satisfied by (|51() thanks to the 

orthogonality of the V>n's, i-e dx ip n (x)ipk(%) — 5k,n- Within simple algebra, that last expression can be rewritten 
as a simple sum [l7j 



^exp 



$(a,a',y) 



v 2 2 



Ma') 



ntp n {o!)ipn(a') + -^A n (a, a') 



nexp 



y 2 2 



(52) 



with 



A n (a,a') = - 



f 7T "I / O/ — Cx 

{I- a - 2a) sin n-(a + a + 2) + 1 + a - 2a - 2- | sin 

- 2 J V \a — a 



n—(a — a 
. 2 V ' 



(53) 



In formula (|48|) . third variable y of the function $ is proportional to Ad/£. Thus the low and high temperature regimes 
correspond to the limits y — > oo and y — > respectively. In the next subsections, we derive asymptotic formulas for 
$ in those two limits. 



2. Low temperature regime 

The asymptotic form of $ when y — > oo limit is obtained by keeping only terms n — 1 in formula (|52|l . This leads to 

*(a,c/,y) = rt{a<) (l + l A ( l{ ^i1\? l +0(ex P [-y 2 /2]) (54) 
V ?T Vi(")'0i(a ')/ 

= cos 2 (?f) [l + ^(a'-l)tan(^') 

(IT (y! \ / 7T \ /" 7T \ ~ 7T 

— J [(M Q , Q - - l)tan^-M Q , Q -J + (1 + m a , a .) tan (^-m a , a >JJ ^ + O (exp [-y 2 /2]) , (55) 

where m a ^ a ' = min(a, a') and M a ^ a i — max(a, a 1 ). Notice that the leading the term is normalized to unity (in other 
words, the leading term satisfies the normalization condition of $). Fig. [S] shows the comparison between asymptotic 
expression <|55|) and the exact formula l">l't evaluated numerically. Asymptotic formula is really accurate for any a, 
even for y close to unity. 
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FIG. 5: Dimensionless function $(0,0?', 2) for a — —0.99 (circles), a — (plus signs) and a — 0.99 (squares). Solid lines follow 
from the numerical calculation of expression 1521 . while points correspond to asymptotic formula 1551 . 



3. High temperature regime 
In order to obtain the small-?/ behaviour of the Poisson transform (|90[1 is applied to 1)52(1 . We find 



$(a,a',y) = 



(1 — a) + (1 — a) exp 

1 1 nM 

2 + 2 



y 



-n(2n — a — 1) 



nexp 



(1 + a) + ^(1 + a + An) exp 

^2 



r n(l + a + 2n) 



2y 



■(a! + a- 2n)(a' - 1 - 2n) 



nexp 



2t/ : 



-(a' -a-2n)(l + a' - 2n) 



(1 + a) + (1 + a + An) exp 

n#0 



-— n(l + a + 2n) 



(56) 



for a' < a, while 



(1 + a) — 2 exp 



${a,a',y) 



1 1 

2 ~ 2' 



nexp 

n#0 



._(a'-a)(l + a') 



(1 + a) ^ exp 



-n(2n — a — 1) 



(1 + a) + ^(1 + a + 4n) exp 



-— n(l + a + 2n) 



7 (a' + a- 2n)(a' - 1 - 2n) 

2y z 



2J(1 - n)exp 



? (a' - a - 2n)(l + a - 2n) 

2y 2 



(1 + a) + (1 + a + 4n) exp 

n#0 



-— n(l + a + 2n) 



^37) 



for a' > a. For y small, numerical estimates of those expressions are obtained by truncating all sums to terms 
n = ±1. They are compared with exact formulas l(56|l - l(57|l in Fig. Eland are quite accurate for any a, and even for y 
close to unity. Let us notice that, in the very high temperature regime and for almost all values of a, $ is close to 
e(a-a')/(l + a). 
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"-1 -0.5 0.5 1 -1 -0.5 0.5 1 

a' a' 



FIG. 6: Dimensionless function <I>(a,a', 1.5) (left panel) and $(a,a',0.5) (right panel) for a = —0.99 (square) , a — (cross) 
and a — 0.5 (circle). The solid lines follow from numerical calculations of formulas 156H - I57|I . while symbols correspond to 
small-y asymptotic formula. 



V. APPLICATIONS OF THE ERGODIC APPROXIMATION FOR CONFINING POTENTIALS 



In the remainder of the paper, we will restrict ourselves to a symmetric confining potential. Let us focus at first on 
the following question: does the ergodic approximation lead to reliable estimates in the low temperature regime ? 
Fully analytical results will be presented for the simplest potentials, namely the harmonic potential and the infinite 
square well. For more complicated potentials, we combine both analytical and numerical methods. 



A. Low temperature regime 



1. General asymptotic formula 



Inserting the asymptotic expressions (|55|l and H89fl for and g (respectively) into formula l|48|) . we find the low- 
temperature behaviour of the ergodic density matrix, 



with 



Perg ^3 P) 



E{i) = 



a) 



-00 2\2 

dll^D sin 2 



8£ 4 



ir(x+£) 



21: 



exp [-(3E(£)] (exp [-<p+ {£, x/l)] + exp [-tp- {I, -x/£)] ) , (58) 



8m£ 2 
2ml 



da' cos 2 )V(a'£), 



(59) 



lm£ 2 f + \ .. , . /7T , 
— =- / da (a — l)tan — a 
7r/r J_i \2 



7ra 
~2~ 



2m£ 2 
irh 2 

2m£ 2 
■kK 2 

2ml 2 



da'V(±a'£)cos 2 ^- 



V(±a'£) 

(M ata > - 1) tan (^M a ^ + (1 + m a , a <) tan y~m a , a <) 



2 / da'V(a'£)sm(TTa')a' - / da'V(a'£)sm(ira') 



■a tan 



/7ra\ 
\~~2~) 



da'V(a'£) cos 2 



7ra 
~~2 



dz V{a'£) cos 2 



(60) 
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Positivity of ip^ enforces the convergence of the integral in the r.h.s. of (|58|) (contributions from boxes with large sizes 
£ do vanish in an integrable way). Expression 158fl can still be simplified in the low temperature regime (/3 — > oo), 
by using the saddle point method. The divergence of potential V(z) for z large ensures the existence of a value £q 
which minimizes the function E{£). We emphasize that £q depends, of course, on the potential, and is attained in 
integral (|58l) only for values of x such that \x\ < £q. Assuming that the second derivative of E is well defined at £o, 
we can apply the saddle point method for \x\ < £q. This provides the general result 

p erg (x,x,f3) (3 ~°°n(x)r(f3) (61) 
for \x\ < £q. In i|61|l . n(x) is the (unnormalized) density 

(exp [-<p+(io,x/£ )] +exp [-y- (£ - x /£ ))) , (62) 

while the temperature dependence is entirely embedded into 

We stress that temperature and position dependences are factorized in formula Ij61|l , like in the exact low-temperature 
behaviour (|17fl . Although the temperature dependence of the ergodic density matrix is not entirely correct (because 
of the presence of factor Xo in Y{j3)), the quantity E(£q) can be identified as the ergodic groundstate energy: indeed, 
it controls the exponential decay of p erg (x 1 x, (3) when [3 — > oo. Below, we show, through several examples, that E(£q) 
is a good approximation for the real groundstate energy. 



n(x) = — sin 2 
£o 



tt(x + go) 
2£ 



For \x\ > £q, the saddle point £ is outside the integration range of JSHJ. In that case, p erg (x, x, (3) could be evaluated 
by expanding the involved integrand for £ close to Of course, for a given /3, if x becomes sufficiently large, 
Perg{%i x, 0) tends to the classical Boltzmann factor (which, by the way, vanishes exponentially fast). 



2. Infinite square potential 



First, let us consider the infinite square well potential 



{A _ ( -Vo if z < \£ P \ 
V[z) ~ \ +00 if z> \£ p \ 



(64) 



which is the simplest symmetric and confining potential. Identity 

-1 



da' cos 2 



ira 



1 



) tan ^Mc, |(1 'j + (1 + m a ,a') tan (J^m^a^ + (a' — 1) tan (^ a ') = (65) 



allows us to recast expression (|58H as 



( o erg (a;,a;,/?)' 3 ~ 00 9(£ p 



a:)— sm 



tt(x + £ p ) 



exp [/3Vq] exp 



-0 



n 2 h 2 



8m£l 



(66) 



Therefore, the ergodic approximation gives the exact groundstate properties of an infinite square well. That exact 
results cannot be retrieved within usual approximations, like the semi-classical one. 



3. Harmonic potential 
For the harmonic potential V(z) — mu> 2 z 2 /2, we find 
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The minimum of that expression is reached for 



to 



n 2 h / 3 



1/2 



2muj V 7T 2 — 6 



(68) 



and it reduces to 



7T 2 - 6 



1/2 



2 



(69) 



Since the second derivative E"(£q) = Amu) 2 (iT 2 — 6)/(37r 2 ) is well defined, it is possible to apply formula 16111 for 
|a;| < lo- Using 



4(tt 2 - 6) 



-24 + 7T 2 1 



we eventually obtain for \x\ < £q, 



I \ z • 2 

n(x) — — am 



ir(x + £ ) 



2£ 



exp 



4(tt 2 - 6) 



■<i x 2 \ /„ x ( irx \ 



24 + . 2 (^-lj ^3 + 7 r-tan V24 



r(i9) = ^expI-W)]. 



(70) 

(71) 
(72) 



The normalized groundstate probability density, ^(x) = n(x) / J dz n(z) is compared in Fig.[7|with the exact Gaussian 
groundstate wavefunction. 




(mui/h) 



FIG. 7: Probability density of the ground state for the harmonic oscillator. The solid line represents the exact solution while 
points correspond to the ergodic approximation 171H . 



4- Anharmonic potentials 

Now, we consider more complicated potentials V(z) = a n z n /n, where n is an even integer. The associated Schrddinger 
equation depends only on the two parameters h 2 /m and a n . Only one typical energy e p = (ft 2 /m) n ^™ +2 ^a 2 /^™ +2 ' 
and one typical length l p — [fi 2 /(m«n)] 1/ ''' l+2 ' can be built in terms of those parameters. Therefore, the energy and 
the spatial extension of the ground state are proportional to those typical energy and length respectively. Again, 
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n 


2 


4 


6 


8 


10 


E{£ ) - E 

En 


13% 


12% 


10% 


7.5% 


6% 



TABLE I: Comparison of the ergodic groundstate energy E(in) derived by using formula ( I74H . to the exact result En (obtained 
by numerically solving Schrodinger equation). Integer n is the exponent of the anharmonic potential 



formula l|tjl|) can be used to compute the ergodic density matrix elements in the low temperature regime and for 
|x| < In. First, the function E{£) reads 



n 2 h 2 a £ n 



(73) 



where R(n) = J . dz cos 2 (7rz/2)z" depends on n. The minimum of above expression, reached for £q 
£ p (ir 2 /4R(n)y 



,l/(n+2) . 

, IS 



E(£ ) = a n £ n R(n) 



1 1 

2 + n 



(74) 



We also find 

i P ± (£ 0l x/£ ) 



2nR{n) 



2nR{n) \£ 



2 f dz sin(7rz)z™ +1 
Jo 



x/l 



dz sin(7T2;)2; ra 



tan 



2£oJ 



dz cos — 
V 2 



X 



dz cos' 



(75) 



which provides n(x). 



Two important comments are in order. The scaling properties of the groundstate energy and extension are indeed re- 
covered within the ergodic approximation. Table[I]presents a comparison with the numerical resolution of Schrodinger 
equation. The agreement, already good for low values of the exponent n, becomes better when n increases. On the 
contrary, the semi-classical approximation completely fails, since it provides a vanishing ground state energy (see 
Appendix C). Comparisons for the groundstate pdf are presented in Fig. El the agreement is impressive. 



B. High temperature regime 

Either the Wigner-Kirkwood expansion or the semi-classical approximation become well-suited when the temperature 
increases. Does the ergodic approximation provides satisfactory results in the high temperature regime? More pre- 
cisely, what part of the Wigner-Kirkwood expansion is accounted for within that approximation? The corresponding 
calculations are carried out in Appendix D. We obtain 



Perg ? p) 



-0V{x) 



<2t:\ 



D 



p dx 2 12 



(76) 



where Aq is a numerical factor close to 0.8. Thus, the leading classical term is indeed recovered, while part of the 
fi 2 -correction is also correctly reproduced. At that order h 2 , the discrepancy of (|76|) with exact expansion i|14[l . 
arises from correlations between occupation times which are not taken into account in the ergodic approximation. 
Nevertheless, that approximation turns out to be also quite reasonable at high temperatures. 



C. Intermediate temperature regime 

In this last section, we present results obtained trough the numerical implementation of formula l|48|l at intermediate 
temperatures. As shown below, the ergodic approximation is particularly efficient in that regime, or in other words, 
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it provides a good interpolation between the exact behaviours at high and low temperatures. Here, we consider 
thermodynamics quantities, i.e the internal energy U = —dp (log (Z(f3)) and the probability density function (PDF) 
^(x, (3) = p(x,x,(3)/Z (0). Comparisons are made with either exact results or other familiar approximations. 

The ergodic approximation appears to be rather accurate for estimating the considered thermodynamic quantities. 
For an harmonic potential, it is well-known that the semi-classical approximation turns out to be exact. As shown 
in figure even in that case, the ergodic quantities are then close to the exact ones. For anharmonic potentials, 
there always exists a large-temperature domain where the ergodic approximation is significantly better than other 
considered approximations. FigureHUIshows results for the potential V(x) — x 1Q . Above temperature domain becomes 
larger when the anharmonicity of the potential increases. 
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FIG. 9: Relative error for the internal energy (left panel) and mean-square error for the probability density function (right 
panel), as functions of the temperature for the harmonic potential. Circles correspond to the ergodic approximation, square 
refers to Wigner-Kirkwood expansion truncated up to order h 2 and diamonds are the classical values. Quantity Ao,i is the 
energy gap between the ground state and the first excited one. 




FIG. 10: Relative error for internal energy (left panel) and mean-square error for the probability density function (right 
panel), as functions of the temperature for the potential V(x) — x 10 . The exact results are obtained by using the spectral 
decomposition ©. Circles correspond to the ergodic approximation, square refers to Wigner-Kirkwood expansion truncated 
up to order h 2 , crosses describe the semi-classical approximation and diamonds are the classical values. Quantity Ao,i is the 
energy gap between the ground state and the first excited one. 
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oo 

E 



sin(x) 



(it — x) 



tan x 

cos (a; 
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, ~ix 6 [-27T, 2tt[. 
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Appendix A: Statistical weights of boxes 

Eigenstates and eigenvalues of the auxiliary Hamiltonian given in formula l|18|) are obtained by solving the 
Schrodinger equation 



d 2 $ _ 2m V 



dz 2 



h 2 



l-Q(z + £) + ®(z-£) 



E 

To 



(79) 



The solutions are obtained by distinguishing three different regions. After introducing dimensionless function ip(z) 
y/l$(z) we obtain 



Aexp[K(z + i)] , if zs^-£ 
ip(x) = { Bsw[k(z + £) + 56], VzG [-£,+£] . 
Cexp[-K(z-£)}, ifz^e 

The continuity of the wave function and its first derivative in ±£ leads to the following four conditions 

k 

A = B sin (59) = —B cos (59) , 
K 

k 

C = B sin [2k£ + 59] = B cos [2k£ + 59] , 

K 



(80) 

(81) 
(82) 



where we have introduced both wavevectors k = (2mE) 1 / 2 /h and K = k(Vo/E — l) 1 / 2 . Rather than deriving the 
general expression for A, C, k, 59 as a function of B, it is sufficient to determine the asymptotic expression for large 
values of potential barrier Vq. In the infinite- Vq limit, the dominant terms are 59 ~ k/K, A ~ Bk/K, k ~ nir/(2£) 
and C ~ (— l)( n ~ x >Bk/K . Finally, the value B = 1 is obtained through the normalization condition. 

Once those results are introduced in (|80|) . if we define the eigenvalues E n = Tr 2 h 2 n 2 / (8m£ 2 ) = E\r? of the n-th 
eigenstates in the large- Vb-limit, we then obtain the dimensionless eigenfunctions of the system in that asymptotic 
case 



Vo—*oo 



-n. 



ifz=-£ 



— (z + £) , Vze] -£,+£[ 



(_l)(»-i) A /r± n , if j i, 

Vq 



(83) 



By using the spectral decomposition JSJ), we obtain the density matrix elements of the free particle constrained to 
stay in the box [—£, +£] 



9° e (-£,x,p) = p° e (x,-£,p) 



p° (£,x,p) = p° e (x,£,0) 



Vb— >oo 1 / E\ v ^ 
V ->oo 1 E- 



n sin 



— (z + £)\ exp [-0E& 2 ], 



n=l 



TITT . 

2£ ' 



which leads to the following expression for the statistical weight functions 

x 7tAd\/s 



g± (x,£,s,f3) 



7T 2 A 3 



t 2£ 



where 



Sf(a, y) = £ n ^ l T^ sin + 1 



•5',; 



exp 



x 7tAdV1 — s 
7' 2l 



2/ 2 a 
— n 
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(84) 
(85) 

(86) 
(87) 



Those two functions satisfy the equality <? + (£, £, s,/3) = g (—x, £, s, 0), the normalization condition 



d£ / ds[g~(x,£,s,f3)+g + (x,£,s,(3)] = 1, 



(88) 
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and they take a simple asymptotic form in the small-temperature regime 



g±(x,£,s,P) X ^ >1 V^^sin 2 
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r tt 2 \ 2 i 


2t 


exp 


8£ 2 



(89) 



The asymptotic expression of g± in the high-temperature regime is obtained by taking the Poisson transform of 
expression 1)87(1. According to Poisson formula, if the Fourier transform of a function / exists, then, for all values of 
Ay 



cx> +oo »+oo 

Ay ]T f(nAy)= £ / dzf( 

^ ™ J OO 



z) e 



-vml'KZ I Ay 



(90) 
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Application of that transformation to expression (|87|) provides 

3/2 



sf(«.i/) 
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Therefore, in the high temperature regime, g± reads 

X D «i [2 I 2 



g (x,£,s,f3) = g + (-x,£,s,P) 



(1 + x/lf 



■ exp 



e {i + x/if 



2X1 s(l-s) 



(92) 



Appendix B: Derivation of the mean occupation time by using the operator method 

Expression (|52|1 for the averaged occupation time can be alternatively computed by introducing the effective 

Hamiltonian H x > — H® + V x >, with (x|V^< = e8(x — x')(x\. Indeed, using a Taylor expansion of the exponential, the 
density matrix elements associated to that hamiltonian can be first written as 



(Xi | exp [-PH X >] \x f ) _ (xi\ exp [~(3H°] \x f ) + e( Xl \d e [exp [-/3H x ,]} e=0 \x f )+0(e 2 



( Xi \ exp [-0H$]\x f ) (xi\ exp [-[3H°]\x f } 

(xi\d e [exp[-PH x ,]] e=0 \x f ) , 2 x 



(93) 



(xi|exp[-/3i/ £ u ]|a;/) 



The Dyson derivation formula l|23|) and the introduction of a closure relation lead to 



ds / dz ( Xi \ exp [-(3(1 - s)H$] \z)S(z - x)(z\ exp [-/3sif°] |x/) 



(x,|exp [-0H x ,]\x f ) _ J0 .,, , 2 . 

(x,| exp [-/3#°]|x,) P£ (^| exp [-/3fl?] |x,) + U [£ ' ^ 

f ds (xi\exp[-p(l-s)H$]\x')(x'\exp[-psH$]\x f ) 
= — f jyr, — ; + 0(e 2 ). (96) 

Above relation can be further simplified by introducing the wave functions <j)%,{x) = £ -1 ' 2 ^ n (x), given in Eq. I|83[l . of 
the constrained Hamiltonian and its associated eigenvalues E n — E\v? . We obtain 

( Xi \exp[-PH x ,]\x f ) _ Pe Jo \ k / \ „ L + (e 2 ) (97) 
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The left-hand-side can be rewritten by using the Feynman-Kac formula as 



(xi\exp[—PH x ,]\xf) _ J 


' £V(£) ex P 


r" £ J 


ds 8(x' — Xi(l — s) — sxj — \d£{s)) 




(xi\exp[-0H°]\x f ) ^ 


UJ 



(99) 



(The effect of part H® of the Hamiltonian is to restrict the integration domain to the subset cu.) Finally, using a 
Taylor expansion of the exponential leads to 



[xj\ exp [-(3H x >]\x;) 
(x t \exp{-pH°]\x f ) 
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= l-pe (9(x'))^+0(e 2 ) 
Combining Eq. J§3J and (|102|l leads to the final result lO . 
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Appendix C: Semi-classical approximation for potentials V(z) = a„z n in the low-temperature regime 

In this appendix, we present the semi-classical analysis for potentials V(z) = a n z n , where n is an even integer and a T , 
a positive real constant. The semi-classical density matrix elements are given by Dashen Formula [l8j 



p sc (x,x,(3) 



exp 



dP (x,Xf) 



dx f 



1/2 



(103) 



where S c {.] is the classical Euclidian action of the trajectory starting from the initial point x and coming back to that 
point in a time (3h, and Po(xi,Xf) is its initial momenta. 



The classical Euclidian action for the particle moving in potential —V 

r /3h 



S r \z\ = 



dt 



jZ 2 + V(z(t)) 



(104) 



can be rewritten, if one introduces the classical energy E — (mz 2 )/2 — V(z) and uses that the exponent n is even, as 



S c [z] = phE + 2 / dtV{z(t)). 



0h 



(105) 



For negative energy values E, we introduce characteristic lengths z m = (—E/a n ) 1 / n and £ p = [h 2 j ' (ma n )] 1 ^ n+2 ^ and 
characteristic energy e p = an£ p . This leads to 



9 . \ 1/2 / „ \ n/2 



1/2 



Expression (|105|l can thus be rewritten as 



S c [z] = I3he p 



-1 + 2, 



/_2_ 

(3e p 



Xd 



I \ n/2 f x x ( ™ +2)/2 



(106) 



(107) 



with F n (y) = y -(«+2)/2) $y&uu n /y/u n - 1. The turning point Zjyi is determined through the constraint that the 
duration of the motion from z m to x is precisely [3h/2. 
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In the more general case, if Xi = z(Q) is the starting point of the trajectory and x/ = z(f3H) is it's ending point, the 
constrain over time is 



(3h = 



Using expression (|1U6|I . that constrain can be rewritten as 
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(108) 



(109) 



with I n (y) = du /\/u n — 1. In the low temperature regime (i.e. when fih tends to infinity), we can derive the 
asymptotic expression of the characteristic length z m . Setting zq = z m (x,x) the solution of implicit equation l|109|) 
for Xi =Xf = x, we calculate z m (x,Xt) — z + 5z if Xf = x + Sx is very close to x. The linearization of equation l|109|) 
then leads to 



Sz Sx 



z o x 2 + (n - 2) (^) I n (x/z ) {(x/z ) n - 1] 



1/2 
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The quadratic case n — 2 is singled out since equation (|109fl can be analytically solve for all values of (3. We obtain 

(111) 



cosh 



and S r \.] = mx 2 \ — 2 - tanh [ —a 



m 



2 V ni 



which ultimately provides the exact density matrix. The non-quadratic cases (n > 2) are more interesting. At low 
temperatures, we find 



z m ^ £p 



l/(n-2) 



and S c tt~2\x\y/2^(^j ' F n 



(112) 



where /„ = lim^^oo I n (y) and F n — lim^oo F n (y) are numbers of order 1, estimated numerically (see table [nj. 



n 


4 


6 


8 


10 


F n 


0.33 


0.25 


0.20 


0.16 


In 


1.31 


0.7 


0.48 


0.37 



TABLE II: Value of F n and I„ for different even exponents n of the potential V(z) = a n z" 



To achieve the calculation of formula (|103fl . we have to compute the partial derivative of the initial momenta with 
respect to the ending point of the classical path. Energy conservation reads 



P (x,x f ) = (2m £p ) ( — 



— 1-1 

z. 



For z m close to zq, we linearize Po(x,Xf) = Pq{x, x) + 6Pq(x, Xf), with the result 
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In the quadratic case, we find 



dP (x,x f ) 



dxf 
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(116) 



For non-quadratic cases (n > 2), we have the low-temperature behaviours 
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Finally, by combining 1103(1 . 1(116(1 and 1(111(1 . we obtain the exact density matrix for the harmonic potential 

(2ma 2 ) 1 /4 exp 



p sc (x,x,/3) 



2 / 2 



^tanh(£V27^i 



(2 7 rfi) 1 /2 S i n h [(/3n/m)V2mai] 

For anharmonic potentials, the semi-classical density matrix behaves as 



(118) 
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(119) 



When /3 — > that expression does not reproduce the expected exponential decay with respect to temperature. 



Appendix D: Ergodic density matrix at high temperatures. 



Let us focus on the calculation of p erg {x,x, 0) in the high-temperature regime. We obtain it's the asymptotical form 
for values of x such that |o;|/A.d 1. For simplification purposes, we restrict our analysis to negative values of x. In 
that case, using asymptotic expression 1(92(1 . we realize that we can omit the paths taken into account by the density 
of measure g + {x, £, s, (3) since their contributions are negligible (on the contrary, for positive values of x, we can omit 
the contribution of paths implied in g~ (x, £, s, /3)).Thus, expression ((48(1 reduces to 
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That expression clearly shows that paths with extension larger than A u with respect to x give a vanishing contribution 
to the integrals. Since we are interested in the limit |a;|/A.D ^> 1, we see that for leading paths x / £ = —1 + Aa with 
Aa <C 1. By starting from 1(56(1 - ((57(1 . we find the asymptotic form of the function $(a, a', y) in that parameter range, 
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(121) 



$ 2 (Aa,Aa',y) = ^ ex P ( ^^ Aa ' 2 ) sinh ( ^2 AaAa ' ' A " A " 
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Then, the integral appearing in the exponential of 1(120(1 can be calculated by using a Taylor expansion of the potential 
around the position x. We obtain 
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V (a'*) = K(a;) + 2 s) + F fc (£, 1 - a)] (122) 



fe=i 



with 



V r fc (^,u) = u 



t k d k V 



jfd** (l + Aa,-l + z,^|^) ( Z -Aa)\ 



(123) 



25 



The explicit forms of the first terms are 
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and V 2 (£,u) = 
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where erfc(z) is the complementary error function. Expression H12Q(I can be formally rewritten as 



p erg (x,x,0) ~ I dt 
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fc=i 



Since the temperature is high, part of the exponential implying derivatives of the potential can be expanded in power 
series. The truncation of these series up to second order with respect to Xd provides 
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After computing the remaining integrals over s and £ we eventually obtain 
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with 
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a numerical factor close 0.8 (estimated numerically to 0.796). 



